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Two Epochs of Globular Cluster Formation from Deep 
Fields Luminosity Functions: Implications for Reionization 
and the Milky Way Satellites 
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ABSTRACT 

The ages of globular clusters in our own Milky Way are known with precision of about 
±1 Gyr, hence their formation history at redshifts z £ 3 and their role in hierarchical 
cosmology and the reionization of the intergalactic medium remain relatively unde- 
termined. Here we analyze the effect of globular cluster formation on the observed 
rest-frame UV luminosity functions (LFs) and UV continuum slopes of high redshift 
galaxies in the Hubble Ultra Deep Fields. We find that the majority of present day 
globular clusters have formed during two distinct epochs: at redshifts z ~ 2 — 3 and 
at redshifts z ^ 6. The birth of proto-GC systems produce the steep, faint-end slopes 
of the galaxy LFs and, because the brightness of proto-GCs fades 5 Myrs after their 
formation, their blue colors are in excellent agreement with observations. 

Our results suggest that: i) the bulk of the old globular cluster population with 
estimated ages £ 12 Gyr (about 50% of the total population), formed in the relatively 
massive dwarf galaxies at redshifts z ^ 6; ii) proto-GC formation was an important 
mode of star formation in those dwarf galaxies, and likely dominated the reionization 
process. Another consequence of this scenario is that some of the most massive Milky 
Way satellites may be faint and yet undiscovered because tidal stripping of a dominant 
GC population precedes significant stripping of the dark matter halos of these satel- 
lites. This scenario may alleviate some remaining tensions between CDM simulations 
and observations. 

Key words: (Galaxy:) globular clusters: general, galaxies: luminosity function, mass 
function, cosmology: theory 



1 INTRODUCTION 

Although globular clusters (GCs) are well studied in our 
local galactic neighborhood, their role in hierarchical cos- 
mology and their formation history at redshifts z > 3, when 
the universe was less than ~ 2 Gyrs old, remains an un- 
solved problem. Several groups have attempted to model 
the formation of GCs using cosmologica l simulations (e.g., 
iGriffen et all l20ld : IZcmoozi et al.1 l201lt) . however, insuffi- 
cient spatial resolution and a poor understanding of their 
initial mass function, as well as the mechanisms responsible 
for their formation leave many unanswered questions about 
their cosmological origin. 

Nearly all GCs are compact self-gravit ating star clus - 
ters fairly homogeneous in heavy elements l|Snedenll2005h . 
This suggests that most of their stars formed in a nearly 
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instantaneo us burst with high efficiency of gas- to-star 
conversion (Uames et al ] 12004 ICarretta et al l 12009 . GCs 
are often classified into two categories: metal-poor with 
[Fe/H]< -1.5 and metal-rich with [Fe/H]^ -1.5. The 
origin of these two populations is unknown, however this 
bimodal prope rty has been observed in multiple galactic 
environments (|Zepf fc Ashman! 1 19931 ). The ages of many 
Milky Way GCs and a few extra-galactic GCs have been 
estimated by analyzing the color-magnitude diagrams of 
the individual GCs. Vertical, AV, methods tend to have 
error bars of ~ 1 Gyr due to uncertainties in defining the 
precise main sequence turnoff, while horizontal methods, 
A(B-V), have similar error bars as a result of uncertainties 
in the theoretical model for the effective temperature 
T e jf of the stars (Saraicdini 200l|). In addition, galactic 
metal poor GCs are older than metal rich GCs and have 
an age spread of approximately 1 Gyr, c ompared to the 
~ 6 Gyr spread for m etal r ic h GCs I Chabover et al.l 
1 19961 : iRosenberg et~ai] 1 19991 : ISalaris fc Weiss! 120021 : 
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|Pe Angeli et all 120051 ; iMarin-Franch et al.l 120091 ). The 

age gap between metal-poor and metal-rich GCs is greater 
than the range present w ithin each individual population 
l|Marin-Franch et al.ll2009T ), suggesting that there are two 
distinct epochs of GC formation, adding to the complexity 
in modeling their formation and evolution history. 

Assuming that the bulk of the old GC population 
formed around the epoch of reionization, iRicottil (j2002T ) 
has shown that they would dominate the reionization pro- 
cess. Ricotti's argument is simple: from the census of GCs 
in the present universe (estimated from the specific fre- 
quency of GCs), he derived the fraction of baryons in GCs 
and the number of ionizing photons per baryon per Hub- 
ble time they emitted (assuming they formed over a 1 Gyr 
period at z ~ 6). This number is likely to be a lower limit 
since most GCs have been destroyed or have been stripped 
of stars due to dynamical processes and stel lar evolution 
l|Ostriker fe Gnedinl 1 19971 ; iFall fc Zhanell200ll ). Today, the 
fraction of stars in GCs today is negligible when compared 
to all of the stars in bulges and disks. However, at red- 
shift z ^ 6, only a small fraction of today's stars existed, 
while many of today's GCs, being the oldest known stel- 
lar systems, may have been present. Ricotti concluded that 
proto-GC formation may have been the dominant mode of 
star formation in primordial dwarf galaxies. In addition, 
because proto-GCs formed with a high star formation ef- 
ficiency and are observed in the outer parts of dark matter 
halos, the UV radiation they radiate can more easily es- 
cape into the intergalactic medium (IGM) l|Ricotti fc Shulll 
I2000T I. and contribute to the reionization process. Indeed, 
assuming a Salpeter IMF and UV escape fraction f esc ~ 1, 
Ricotti showed that GC formation at z ^ 6 can easily pro- 
duce enough UV photons to reionize the universe. 

The main uncertainty in this model is the fraction of 
the observed population of old GCs that formed at red- 
shifts z ^ 6. The epoch of reioniza tion likely extended from 
z ~ 14 to 2 ~ 7 (|Fan et al.ll2006h . The uncertainty on the 
age determination of old GCs prevents us from discriminat- 
ing whether old GCs formed at redshift z ~ 3 — 4, after 
the epoch or reionization, or at redshift z ~ 7 — 14, thus 
contributing to reionization. In order to understand the role 
that GCs may have played during the reionization epoch, 
the formation history of the old GC population with esti- 
mated ages ^ 12 ± 1 Gyr (i.e., formation redshift z £ 4) 
must be better constrained. In addition, constraining proto- 
GC formation at high redshift may be a key element to un- 
derstanding the o rigin of the classical and ultra-faint Lo- 
cal Gr ou p dwarfs feicotti fc Gnedinl [20051 ; iBovill fc Ricottil 
l2011aH bl; [Bovlan-Kolchin et alJl201ll ). 

Figure [T] (left) shows the distribution of the redshift of 
formation of the Milky Way GCs b a sed o n age measure- 
ment compiled by iForbes fc Bridges! |2010). The distribu- 
tion shows that about 55% of Milky Way GCs are older than 
12 Gyrs, thus formed at redshift z > 4. However, the errors 
on these ages are on the order of ±1 Gyr. Assuming the 
worst case scenario where all of the GCs are 1 Gyr younger 
than the data suggests, 20% of Milky Way GCs would have 
still formed after redshift z > 4. The right panel of Figure [T] 
shows a Monte Carlo realization, consistent with observa- 
tions of the redshift distribution (including Gaussian errors 
on the ages) of a bimodal population of GCs that form 20% 
of all clusters at 3 < z < 3.5 and 80% at 5 < z < 5.5. Due to 



the large errors on the age measurements, is impossible to 
discriminate between different choices of the underlying GC 
bimodal distribution. For instance, a bimodal distribution 
with 45% of GC forming between 2.5 < z < 3 and 55% at 
6.5 < z < 7, is statistically indistinguishable from the one 
shown in the figure. 

The old stellar populations and compact size of present 
day GCs adds to the difficulty in detecting them beyond the 
Milky Way, even w ith our most powerful space telescopes. 
iKalirai et al.l (|2008l ) has identified bright unresolved sources 
surrounding the parent galaxy at z = 0.0894 in a Hubble 
space telescope (HST) image, which is one of the furthest 
known systems of old GCs at a luminosity distance of about 
404 Mpc. 

For about 5 — 10 Myrs after their formation, individual 
proto-GCs are extremely bright and have blue broad-band 
colors. In principle they are detectable with current instru- 
mentation up to redshift z ~ 7 (see Section |3J) and are cer- 
tainly within the reach of even deeper near-IR observations 
which wil l be achieved with the James W ebb space telescope 
(JWST) (|Schaerer fe CharbonnefeoilT ). Moreover, GCs are 
typically observed as sy stems of 3 — 10 GCs in dwarf galaxies 
(Puzia fc Sharinall2008l ). to hundreds in more massive galax- 
ies (the Milky Way has about 150 GCs and Andromeda 
about 500). With current space telescopes, proto-GC sys- 
tems forming in high redshift dwarf galaxies are poorly re- 
solved spatially. In Section l3.2l we show that they appear as 
blurred extended objects, brighter than an individual GC, 
since their formation is likely nearly coeval given the typical 
dynamical time scale of a few tens of Myrs in high-z dwarfs. 
These young proto-GC systems are therefore more easily de- 
tectable than sparsely distributed individual GCs in larger 
galaxies. 

In this paper, we take an indirect approach to identify- 
ing the redshift of formation of GCs. We seek to constrain 
the largest number of proto-GCs that are allowed to form 
within any given redshift range, while remaining consistent 
with the observed UV luminosity functions (LFs) and col- 
ors of galaxies in the Hubble Ultra Deep Fields. We also 
produce simulated images of young proto-GC systems to 
lay the foundation for their direct identification with Wide 
Field Camera 3 (WFC3) on the HST and the Near Infrared 
Camera on the JWST. 

This paper is organized as follows. In Section [2] we 
present basic calculations of the luminosity and mass bud- 
get of proto-GCs. In Section [3] we discuss the capabilities 
of the HST and the JWST in detecting and resolving young 
proto-GCs and present synthetic images of proto-GC sys- 
tems based on o bservations of old GC s ystems in nearby 
dwarf galaxies by IPuzia fc Sharinal |2008T ) . In Section [4] we 
constrain the proto-GC formation rate at high-z by compar- 
ing observed LFs to simulated LFs of forming proto-GCs, 
and calculating the effect of proto-GCs formation on the 
UV continuum slope at 1500 A. We present our results and 
discussion in Section [5] and the summary and conclusions in 
Section [6] . 
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Figure 1. Left. Distribution of the redshift of formation of the Milky Way GCs based on age measurement from Forbes & Bridges (2010). 
The distribution is poorly known at z > 3 because the typical error on the ages of GCs is ±0.5 — 1 Gyr. Right. Monte Carlo realization 
of the redshift distribution (including errors on the ages) of a bimodal population of GCs that form 40% of all clusters at 3 < z < 3.5 
and 60% at 5 < z < 5.5. Due to the large errors on the age measurements, it is impossible to discriminate between different choices for 
the underlying GC bimodal distribution. Most choices, as the one shown here, produce observed distributions broadly consistent with 
the Milky Way population of GCs. 



2 BASIC CALCULATIONS AND 
ASSUMPTIONS 

The mass distribution of present day GCs in the Milky Way 
is well described by a logno rmal distributio n with typical 
mass Aiobs ~ 1.4 x 10 5 M© (|Shu et alj|2010t ). However, the 
initial masses of young proto-GCs, Mint, were significantly 
larger than today's observed masses, M obs . Stellar evolution 
and dynamical effects, including galactic tides, dynamical 
friction and two-bo dy relaxation dOstriker fc Gnedinl 1 19971 ; 
iFall fc Zhang] I 200 ll ). significantly reduce the number and 
mean mass of young GCs from the epoch of formation to the 
present. Assuming a Sa lpeter initial mass function (IMF), 
IPrieto fc Gnedinl (|2008l ) found / di = M ini /M obs ~ 9.1. In 
their simulations they include the effects of stellar evolution, 
two-body rel axation, dynam i cal fr iction and tidal shocks as 
described by IFall fc Zhanel ((2001). Simulations also show 
that the current mass distribution of GCs can be repro- 
duced if the GC initial mass function (GCIMF) is either a 
power-law dN/dlogM oc Mr° A7 with 10 4 < M ini < 10 7 M 
or a lognormal distribution wit h characteristic mass Mini = 
1.5 x 10 6 M Q (|Shu et al.ll20l"oh . 

Although the destruction rate and mass loss of GCs 
remains somewhat uncertain because these quantities are 
based on theoretical models with poorly known initial con- 
ditions (such as the IMF of stars in GCs, the GCIMF, 
and the formation sites of proto-GCs), recent observa- 
tional developments on the stellar populations in GCs pro- 
vide independent evidence fo r values of / d i ~ 10 — 25 
l|Schaerer fc Charbonnell 1201 ll ) . It is now clear that nearly 
all GCs exhibit large star-to-star variations of abundances 
in light elements (from C to Al), interpreted as a signature 



of some degree of self-enrichment by gas pr e-processed in 
short -l ived, h igh mass stars (for d e tails see. |Gratton et alj 



20011, 120041; iPrantzos et al.l 120071 ; ICarretta et all l2009bl; 
Charbonnel I2010L and references therein). Self-enrichment 



models of GCs assume at least two populations of stars, with 
the first generation of massive stars polluting the second gen- 
eration. The polluters are thought t o be either massive AGB 
stars or fast rotating massive stars (|Prantzos fc Charbonnell 
2006). These models, after assuming an IMF for the first 
generation stars, are constrained to reproduce the present 
proportion of first to second generation stars and their ob- 
served abundance patters. Generally, all the models require 
the initial stellar masses of GCs to be considerably larger 
than their present day values (|p~rantzos fc Charbonnelll2006l; 
Decressin et al.ll2007i iD'Ercole et al.ll200Sl; IDecressin et al l 
20101 : ICarretta et al.H2010h. ISchaerer fc Charbonnell l|201ll ). 



using the model of IDecressin et all ( 20071 ) . have determined 
/di as a function of the IMF of the first generation of stars. 
They found /di ~ 8 — 10 in models in which the second gen- 
eration stars are all retained by the GCs, and /di ~ 15 — 25 in 
models in which second generation stars found in the Milky 
Way halo originate from the present-day population of GCs. 
This result is independent of, but in goo d agreement with 
the e xpectation from dynamical models (|Prieto fc Gnedinl 
2008). Thus, for the rest of this paper we will adopt a con- 
servative fid ucial value /di = 9.1, a t the l ower end of the 
estimate by ISchaerer fc Charbonnel (|201ll ). Note however 
that our results on constraining the formation rate of proto- 
GC systems do not depend on the assumed fiducial value for 
/di. 
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Figure 2. Left. The solid black line represents the average amount of GCs that the NIR camera on JWST would see as a function of 
redshift. The dashed line represents the average amount of GCs that the WFC on HST would see as a function of redshift. This assumes 
a limiting magnitude of 31 for JWST and 30.5 for HST Right. The star formation rate per square parsec as observed with JWST is 
shown as the solid black line and the dashed line represents the same quantity as observed with HST. 



2.1 Luminosity Function of proto-GCs 

For a lognormal GCIMF, with a typical mass of Mini = 
1.4 x 10 6 Mq, a young proto-GC is expected to have a 
peak UV magnitude at 1500 A of Mab ~ —16.7 at an age 
< 5 — 10 Myr. These GCs then fade by mo re than 4 mags 
due to stellar evolution (see e.g. models of iLeitherer et all 
HHl)). For a power la w GCIMF the most massive GCs 
have MZT ~ 10 7 M (|Shu et al.ll20"l0n , corresponding to 
Mab ~ —19.3. At redshifts of z ~ 7 — 10 this would corre- 
spond to a typical UV rest frame magnitude of rriAB ~ 30.3 
at the peak of brightness, slightly fainter than the current 
detection limits of the deepest near-I R images taken with t he 
WFC3 camera onboard of the HST (|Bouwens et al.ll2010T ). 

2.2 Globular Cluster Formation Rate and Mass 
Budget 

Using published data on the spec ific frequency of GCs in 
spiral, elliptical and dwarf galaxies. iRicottil (120021 ) estimated 
that about oj gc ~ fdi(2.7 X 10~ 4 ) = 0.25% of cosmic baryons 
have been converted into GC stars. Although this fraction is 
about 40 times smaller than today's baryon fraction in stars 
(about 10%), it is likely a significant fraction of the col- 
lapsed baryons at redshifts z £ 6. Based on a compilation 
of observations of the p resent number of GCs in galaxies, 
iPuzia fc Sharinal (2008) estimate an efficiency of GC for- 
mation in dark matter halos, r\dm ~ 5.5 x 10~ 5 , defined as 
Mg°* = r/dmMdm (assuming standard cosmological parame- 
ters Qb = 0.04 and S!d m = 0.21). Since at z = most dark 
matte r is in collapsed halos, the results of IPuzia fc Sharinal 
(|2008h imply 0. b /Q. dm ~ 2.9 x 10 4 fdi, which is 

very g ood agreement with the value estimated by IRicottil 
(2002). Hereafter we will adopt a fiducial value for the mass 



density of proto-GCs at formation: 

p gc = uj gc pt ~ 1.35 x 10 7 M Q Mpc" 3 , (1) 

where we have assumed /<& = 9.1 and th e mean cosmic 
baryo n density pi) = 5.51 x 10 9 MqMpc -3 (iKomatsu et all 

Eon]). 

It is easy to show that if even a small fraction f gc ~ 20% 
of today's GCs formed at redshift z > 7 (over a time interval 
At = 0.5 Gyr), the GC formation rate per comoving volume 
would be 

P 9 c « hc^ « 0.54 x 10- 2 AW 1 Mpc" 3 , (2) 

comparable to the star formation rate density observed at 
redshift z ~ 8 (jBouwens et al 1l201ll ). These simple estimates 
suggest that GC formation can be an important mode of star 
formation at high-z, and perhaps dominate the reionization 
process. 



3 DETECTING PROTO-GLOBULAR 

CLUSTER SYSTEMS IN HUBBLE DEEP 
FIELDS 

The planned launch of the JWST provides exciting 
prospects for our understanding of the high redshift uni- 
verse. Here, assuming negligible dust extinction, we address 
the capabilities of both the HST and the JWST to detect 
young proto-GCs at high redshift. First, we estimate the ab- 
solute luminosity, number and morphology of young proto- 
GC systems in HST deep fields and make predictions for 
JWST deep fields. 

The left panel in Figure [2] shows the number of GCs, 
N gc , in the field of view of the Wide Field Camera (WFC) 
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on the HST (dashed line) and the Near Infrared Camera 
(NIR) on the JWST (solid line). For illustration purposes, 
we have assumed that half of all proto-GCs (after correcting 
for destruction), f gc = 50%, have formed in a single episode 
at redshift z. The value of N gc is estimated as 

N gc (z) = 30 A n gc f Um (z)¥^, (3) 

where 9 a is the field of view of the camera, Dl(z) is the 
luminosity distance at the redshift of interest, fu m (z) is the 
fraction of proto-GCs brighter than the limiting magnitude 
of the telescope at that redshift, and n gc is the number den- 
sity of GCs (rigc ~ 39.1 Mpc -3 for a power law IMF and 
n gc ss 33.8 Mpc' 3 for a lognormal IMF). Values of N gc 
for different f gc can be easily inferred by scaling the plot 
in Figure [2] (left). In principle, the WFC3 on the HST has 
the capabilities to detect proto-GCs with Mini = 10 7 Mq at 
redshifts up to z = 7 (see Section f3 . 2 [I . However, such high 
mass clusters are only j mag ~ 0.8% of the total l|Shu et al.l 
2010). N ga peaks at z ~ 7 — 8 because although the num- 
ber of proto-GCs increases with redshift due to the larger 
volume in the field of view, the limiting magnitude of the in- 
strument reduces the fraction of GCs that are detectable to 
the few most massive ones. In conclusion, single proto-GCs 
during their youth are detectable with current instrumenta- 
tion up to z ~ 6, and are certainly well within the reach of 
deeper near-IR observations which will be achieved with the 
JWST. Although we cannot spatially resolve these sources 
at high redshifts, there presence in WFC and NIR images 
may be significant when compared to the number of high- 
redshift galaxy candidates. 

3.1 GCs Versus Host Galaxy Surface Brightness 

Using the average half li ght radius calculate d from the For- 
nax Cluster Survey VII ^Masters et alJl20ld ). of 2.8 ±0.3 pc 
for red GCs and 3.4 ± 0.3 pc for blue GCs, the maximum 
redshifts that WFC and NIRcam can resolve these systems 
are z = 0.0028 (11.9 Mpc) and z = 0.0072 (30.7 Mpc), re- 
spectively. Without adaptive optics from the ground it is 
impossible to resolve or take a spectrum of GCs at inter- 
mediate and high redshifts. Thus, because of their high star 
formation rate at formation, GCs will appear as especially 
bright point sources surrounding their host galaxies and for 
about 5-10 Myrs from formation they may outshine their 
host galaxies. A GC of 10 7 Mq, with a half light radius of 
3.4 pc that forms all of its stars in ~ 10 Myr, will have a 
SFR per unit area of E 9C « 2.75 x 10~ 2 M yr _1 pc~ 2 . How- 
ever, since the proto-GC is unresolved the effective SFR per 
unit area depends on the pixel size of the detector, and at 
redshift z = 6 is ±fj « 1.81 x 10" 5 Moyr^pcT 2 for the 
WFC on the HST. As a comparison, the bulge of the Milky 
Way, which has a mass of approximately 2 x 10 10 Mq within 
a radius ~ 3 kpc, assuming its stars forms over 1 Gyr, would 
have a SFR per unit area of E» « 2.83 x 10" 6 Meyr^pc" 2 . 
The Large Magellanic Cloud (which we are using to com- 
pare GCs to dwarf galaxies) has a mass of 3 x 10 10 Mq 
with a effective radius of 5 kpc. Assuming a continuous 
SFR over 1 Gyr, we infer a SFR per unit area of E* w 
3.82 x 10" 7 M Q yr _1 pc" 2 . Figure fj (right) shows the effec- 
tive SFR per unit area of forming proto-GCs as a function 
of redshift and for two different proto-GC masses. Because 



proto-GCs are unresolved, with increasing redshift their ef- 
fective surface brightness increases due to the increase of the 
angular resolution. Assuming that proto-GCs systems are 
poorly spatially resolved and that their formation is syn- 
chronized to within 10 Myr as described later, we expect 
that the light from a high-z dwarf galaxy will be dominated 
by the proto-GC system for about 10 Myr after their forma- 
tion. 



3.2 Simulated Images 

The best constraint on proto-GC formation at high redshift 
is their direct detection. As discussed in Section [3] the large 
star formation rates and compact sizes of proto-GCs allow 
them to outshine their host galaxies for a brief time. Figure[3] 
shows simulated images of proto-GC systems for 1 Msec ob- 
servations with HST and JWST at different redshifts, and 
for two different choices of the host galaxy mass. To generate 
the images, we have used data on old GC systems observed 
in a survey of nearby sphe roidal and irregular dwarf galax- 
ies (jPuzia &: Sharinall2008l ). We randomly generate systems 
of GCs around these types of dwarf galaxies, modeled af- 
ter the number of GCs, and the magnitude and structural 
data of the host galaxy. We randomly select a radius, de- 
termined by averaging the volumes of shells in which GCs 
orbit around eight different dwarf irregular galaxies. We then 
project on the sky the GC position. The magnitude of each 
GC is also randomly assigned in a range of three magni- 
tude bins from — 18 < Mab < —15. We use the theoreti- 
cal point spread function to convolve the point source over 
many pixels. Next, the three-dimensional positions of these 
systems were projected onto a two dimensional image with 
an arbitrary viewpoint and applied a normalized gray scale 
that attributes relative brightness to the proto-GCs. We in- 
clude the background noise with signal to noise ratios of the 
1 Msec exposure on each camera (consistent with that of the 
HUDF). 

The brightness of the most massive proto-GCs is at the 
detection limit for the two space telescopes. Because of this, 
the signal to noise ratio of these systems is low and it might 
often be impossible to extract these systems over the back- 
ground noise. Additionally, the theoretical images demon- 
strate that at around z = 1, the systems of GCs tend to 
produce an extended source rather than be resolved as dis- 
tinct objects. This may effect the observers ability to dif- 
ferentiate between GC systems and the host galaxy all the 
way up to z — 4. Furthermore, if the systems form in groups 
rather than individually, this effect may be amplified and 
further limit the ability of an observer to identify GC sys- 
tems. However, at z > 4, this effect will be significantly 
reduced and the probability of observing individual GCs or 
GCs in groups is far greater. 



4 CONSTRAINING THE PROTO-GC 

FORMATION RATE AT HIGH-Z WITH UV 
LUMINOSITY FUNCTIONS AND 
CONTINUUM SLOPES 

In this section we calculate the effect of the formation of 
proto-GCs on the observed high-z UV luminosity functions 
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Figure 3. The top rows are the theoretical images from the WFC on HST and the bottom rows are the theoretical images from the 
NIRcam on JWST. These were computed to model a total exposure time of 1,000,000 seconds to be consistent with the HUDF. HST 
sensitivity extends from 2 = — 8 from the UV and IR imaging equipment, while JWST only has sensitivity higher than z = 2.33 for 
1500 A. At the higher redshift we can see that the individual GCs are more identifiable. PSFs for JWST were found in the simulated 
PSF library provided by STScI while the HST PSFs were generated with the Tiny Tim software. Shown at the top left is an actual 
image from Bouwens et al. (2011). We use this to show a comparison of an actual 2 = 8 candidate that appears to be multiple sources 
with the multiple source images that we have generated. 



(LF). The basic idea is to build synthetic luminosity func- 
tions, assuming that the GC system outshines the other stel- 
lar populations in the galaxy during the first 5-10 Myr af - 
ter its formation. We use d ata from [B ouwens ct al.l (|201 lh 
for LFs at redshifts z > 7, I Bouwens et al.| lj2007t) for the 
B, V, and I dropout LFs, and lOesch et al.l (|20T0T ) for the 
z = 1.5, 1.9, and 2.5 LFs. 



4.1 Unresolved proto-GCs 

Our first approach is to consider the formation of individ- 
ual proto-GCs. Although not the most realistic, this case 
is the simplest and nearly model independent. Physically, 
this approach corresponds to a scenario in which the forma- 
tion of the proto-GC system in a dark matter halo is not 
synchronized to within a few tens of Myrs, or the most lu- 
minous proto-GC outshines all the others. This is the most 
conservative assumption but is less appropriate at redshifts 
2 ^ 6 when the typical dynamical times in virialized halos 
are short: 



Rs 



11 Myr 



1 ~f" Zyii 

TO - 



(4) 



where R s — Rvir/5 and the scale radius and the 

virial velocity of the halo, respectively. We have adopted two 
choices of the initial mass functions of proto-GCs (GCIMF): 
a lognormal and power-law (see §[2]), both repro ducing the 
prese nt day GC mass function in the Milky Way |Shu et al.l 
120101 1. Both choices produce nearly identical constraints. 

In order to calculate the proto-GCs LFs, as a fist or- 
der approximation we assume that GCs form in a short 
burst with At form, ~ 10 Myr (instantaneous mode of star 



formation). Using Starburst-99 (|Leitherer et al.l 1 19991 ) we 
derive the evolution of the absolute magnitude of proto- 
GCs and their color (parameterized by /3, the log-slope of 
the spectrum at 1500 A) as a function of time from the 
burst. Proto-GCs produce the majority of their radiation 
in the UV band in the first 5 — 10 Myr of their life, be- 
fore rapidly evolving to higher magnitudes. We compute the 
fraction p gc (M gc ) of GCs that form in each log mass bin, 
using constant log bin spacing of 0.2 Mq, spanning from 
10 4 Mq to 10 7 ' 2 M Q for the power-law GCIMF and from 
10 3 ' 2 Mq to 10 7 ' 4 M Q for the lognormal GCIMF. The num- 
ber of GCs at peak luminosity per unit comoving volume 
is n gc (M gc ) = f g cfonP g cP g c/M gc , where f gc (z), the fraction 
of proto-GCs forming in each redshift bin, is the free pa- 
rameter we wish to constrain. Here, fon is the fraction of 
time the proto-GC is at peak luminosity which is about 
/ on = 10 Myr/5T Z , where ST Z is the time interval cor- 
responding to the characteristic redshift depth of the ob- 
served luminosity function. Finally, we compute the lumi- 
nosity function 4> gc from the proto-GC mass function con- 
verting from stellar masses to magnitudes (using Starburst- 
99 for an instantaneous burst). 

In order to account for the evolution of GCs to higher 
magnitudes as they age, we apply a correction to the calcula- 
tion of the LFs. Proto-GCs after the first 5-10 Myrs become 
fainter and redder, but their number in each redshift bin is 
larger because they remain at fainter magnitudes a longer 
time. We use Starburst-99 to determine the time At a proto- 
GCs of constant mass has a given magnitude. The number of 
GCs at fainter magnitudes from the peak is obtained anal- 
ogously to the calculation above but with f on — At/AT Z . 
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Thus, after converting proto-GCs masses to magnitudes, we 
integrate over a time interval AT Z . 

Figure [4] (left) shows the UV LFs of (individual) proto- 
GCs (dotted lines) compared to the observed LFs in in the 
HUDF at different redshift intervals. The normalization of 
the proto-GCs LF in each redshift bin is proportional to the 
parameter f gc , that we want to constrain. The proto-GCs 
LFs have steeper faint-end slopes than the observed LFs, 
hence we are able to set upper limits on f gc . The values of 
fgc for each redshift bin are reported in Table [S] 



4-1.1 Effective UV Continuum Slopes 

Once we have constrained f gc using the observed LFs, we 
analyze the effect of GC formation on the UV continuum 
slope at 1500 A. Proto-GCs wit h ages < 5 — 10 Myr have 
UV continuum slope /3 = -2.69 jLeitherer et al.lll999T l. We 
evolve the higher mass GCs into the lower luminosity bins, as 
was described when comparing the luminosity functions, and 
compute the weighted average of the f3 at half magnitude 
steps with respect to an arbitrary initial absolute magnitude. 
To determine the effective color, we set /3 e // = f(Pgc—/3)+/3, 
/3g C is the color parameter from the GCs and j3 are the 
observed mean value s of deep field galaxies described by 
iBouwens et~aH l|2009l ). and / is the fractional influence of 
the GCs on the UV continuum slope. In order to constrain 
the maximum percentage of GC formation, we first measure 
/, the ratio of the proto-GC LF <f> gc to the observed LF (f>g a i- 
We then apply this ratio to the corresponding bin in the UV 
continuum slope and determine whether this maximum per- 
centage can be further constrained by a new fraction fp. 
Finally, we multiply the old GC percentage by fp to obtain 
our new constraint. In Figure U (right), /3 e // is compared 
to the the current values of the UV continuum slope. Be- 
cause of the limited range of the UV continuum slope at 
the faint end we extrapolate the least squares fit yet remain 
conservative in constraining fp estimates. When comparing 
/3 e f f to the known UV continuum slope with the Power Law 
and Gaussian models where GCs form in groups of ten, we 
can constrain the maximum percentage of GC formation to 
a higher degree of accuracy by noting where f3 e ff de viates 
from the error bars recorded bv lBouwens et all (2009). 



4.2 Unresolved Proto-GC Systems 

The second approach considers that more than one proto- 
GC forms in each host galaxy, and as discussed in Sec- 
tion [3721 observations measure the integrated light from the 
proto-GC system. In this second case, we do not need to 
know the GCIMF, but the efficiency of proto-GC formation 
as a functio n of the mass of th e host galaxy. We build on the 
results from iPuzia fc Sharinal (2008) that the mass in GCs 
M gCi tot in a large sample of local dwarf galaxies is propor- 
tional to the total host halo mass M dm within the luminous 
radius of the galaxy: M gc ,tot(z = 0) = r/ dm (z = 0)M dm (z = 
0) with r\dm(z = 0) = 5.5 x 10 _J . However, this relationship 
is only measured at z = and does not take into account 
evolutionary effects such as galaxy mergers, proto-GC for- 
mation and disruption as a function of time. For instance, 
we know that only a fraction fdi(Mdm, z) ^ 1 of proto-GCs 



have survived to the present. Thus, we assume that 

M gc ,tot(z) = -q dm (z)M dm (z), (5) 

where n drn (z) is the parameter we wish to constrain. 

Modeling GC systems in dark matter halos is some- 
what model dependent as the shape of the luminosity func- 
tion depends on the range of halo masses in which proto- 
GCs form. However, we find that the only parameter that 
significantly affects the luminosity function is the mass of 
the smallest dark matter halos allowed to form GCs. Our 
model assumptions are as follows. We assume a constant 
efficiency of GC formation as a function of the halo mass 
ri dm (z) for M d ~" < M dm < M"%£ x , and no globular cluster 
formation outside this halo mass range. We use the Press- 
Schechter formalism to obtain the mass function of dark 
matter halos as a function of redshift dnhaio/dM dm , then 
convert dark halo masses to the masses of proto-GCs sys- 
tems using Equation ([5]). Finally, we calculate the luminosity 
function by converting the total proto-GC mass to luminos- 
ity using Starburst 99 as described in § 14.11 As discussed in 
the previous here we also include the effect of aged GCs that 
form within the redshift bin under consideration and in the 
higher redshift bins (note that at high-z there is little time 
for GCs to age significantly and become redder). 

The proto-GC LF is not sensitive to M%£ x , that we 
hence assume larger than the exponential cut-off mass-scale 
in the Press-Schechter mass function. We instead find that 
the results are sensitive to , or equivalently to the virial 

temperature of the halo T™ n . Hence, we explore the results 
assuming 4 different values of 7™" = 5 x 10 4 , 8 x 10 4 , 10 5 K 
and 1.5 x 10 K. In each model we vary the value of r) dm {z) 
to be consistent with observations of the LFs and colors in 
each redshift bin, similarly to what was done for single proto- 
GCs. The main different here is that varying rj dm shifts the 
proto-GCs LF horizontally along the magnitude axis, while 
for single proto-GCs the free parameter f gc had an effect on 
the overall normalization of the LF (shift along the y-axis). 

Figure [5] (left) shows the LF for proto-GCs systems 
(dotted line) compared to observed deep field galaxies LFs 
(solid lines) in different redshift bins. In each panel, cor- 
responding to a different redshift, we have used a value of 
Vdm{z) that maximizes the contribution of the proto-GC lu- 
minosity function to the observed LF remaining consistent 
with observations. We also repeat the same procedure de- 
scribed in Section r4.1.1l to further constrain rj dm (z) from the 
UV colors of observed deep field galaxies. Figure [5] (right) 
shows the UV colors of proto-GCs systems, f3 gc (dotted lines) 
and the effective galaxies UV colors, /? e // (dashed lines), 
compared to observations (solid lines and points with error- 
bars). The curves in the figure are for the largest values of 
Vdm(z) consistent with observations. 

The values of r\ dm {z) constrained by the data allows us 
to calculate p gc {z) = p b rj d m^haios, where ujhaios is the frac- 
tion of collapsed dark matter in halos with M^ n < M dm < 
M"X^ X for each observed redshift bin of depth At(z). From 
p gc we derive the upper limits on the proto-GC formation 
rate p gc {z) ~ p gc (z)/ At(z) and hence upper limits on f gc (z). 
The results are shown in Table [5] for our four model assump- 
tions and are summarized in Figure [6] (left). 
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UV Continuum Slope 
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Figure 4. Left. The UV Luminosity Functions from Oesch et al. (2010) (z = 1.5, 1.9,2.5), Bouwens et al. (2007) (z = 3.8, 5.0, 5.9), and 
Bouwcns ct al. (2011) (z = 6.8, 8) are shown as the thick black line. The observed luminosity function is only shown to the limiting 
magnitude of the observational data. The dotted lines show the LF of a fraction f gc of proto-GCs forming in the same redshift bin as the 
observations, assuming a power law GCIMF. The long-dashed line represents the sum of the proto-GCs LFs and observed LFs. Right. 
The UV continuum slopes for the individual GC formation model. The observational UV continuum slopes are shown with the solid lines 
as the least squares slopes. The least squares fits are only shown with error bars if they were provided by the literature and extrapolations 
are shown with a gap. Raw data points and accompanying error bars are also shown from Bouwens et al. (2009) and Bouwens et al. 
(2011). The dotted lines are the theoretical UV continuum slopes for forming GCs and the dashed lines are the /3 e // slopes. 




Figure 5. The left and the right panels are analogous to the ones in Figure [4] but for systems of proto-GCs allowed to form in halos 
with virial temperature T v i r > 5 X 10 4 K. The long dashed lines show the difference between the observed LFs and the proto-GCs LFs. 



© 0000 RAS, MNRAS 000, 000-000 



Two Epochs of Globular Cluster Formation 9 



Table 1. Constraints on globular clusters formation efficiency from LFs and UV colors. For proto-GCs 
systems we have assumed T v i r > 5 X 10 4 K. 



(*> 


t[Gyr] 


dt[Gyr] 


fgc 


Single GCs 


»7dm(«)/5-5 x 10" 


~ 5 Systems of GCs 


LF 


LF+color 


LF 


LF+color 


1.5 


9.320 


0.99 


75% 




0.23 




1.9 


10.155 


0.92 


100% 




0.79 




2.5 


10.999 


0.56 


100% 


100% 


2.56 


0.77 


3.8 


11.995 


0.36 


55% 


30% 


1.00 


1.00 


5.0 


12.469 


0.28 


35% 


21% 


2.08 


1.04 


5.9 


12.695 


0.17 


40% 


24% 


4.08 


4.08 


6.8 


12.858 


0.13 


28% 


11% 


7.38 


7.38 


8.0 


13.014 


0.14 


15% 




14.46 





Table 2. Upper limits on globular clusters star formation rate (Mq / yr / M pc 3 ) 



z 


T vir > 1.5 x lO 5 ^ 


T vir > 1.0 x 10 5 K 


T vir > 8.0 x 10 4 K 


T vir > 5.0 x 10 4 K 


4.25 


0.0063 


0.0061 


0.0055 


0.0051 


4.75 


0.0048 


0.0050 


0.0035 


0.0032 


5.25 


0.0042 


0.0067 


0.0052 


0.0073 


5.75 


0.0050 


0.0100 


0.0095 


0.0151 


6.25 


0.0061 


0.0105 


0.0118 


0.0181 


6.75 


0.0069 


0.0082 


0.0112 


0.0162 


7.25 


0.0067 


0.0060 


0.0095 


0.0144 


7.75 


0.0055 


0.0051 


0.0075 


0.0146 


8.25 


0.0041 


0.0046 


0.0054 


0.0155 


8.75 


0.0026 


0.0041 


0.0033 


0.0164 



5 RESULTS AND DISCUSSION 

The upper limits on the proto-GC formation rate from the 
Hubble deep fields are summarized in Table [5] and in Fig- 
ure [6] The table shows the constrains from the observed LFs 
and UV colors at 1500 A (parameter 0) of high-z galaxies, 
for both single proto-GCs and systems of proto-GCs. 

Strong constraints are possible because the proto-GC 
LFs have a faint-end slope generally steeper than observed, 
and their spectral slope at 1500 A is typically bluer than ob- 
served (/? ~ —2.7). This is because the instantaneous burst 
of star formation of proto-GCs and the subsequent rapid 
time evolution of their luminosity from bright to faint, steep- 
ens the LF with respect to the intrinsic slope of the GCIMF 
(for single proto-GCs). However, for systems of proto-GCs 
the faint-end slope of the LF depends on the slope of the 
halo mass function, and thus depends on the masses of the 
proto-GC host halos and their redshifts of virialization. This 
is why, for GC systems, the limits are sensitive to the as- 
sumed minimum mass of halos that can host proto-GCs. 

At high-z the limits on the formation of rate of proto- 
GC systems are the most stringent. These, in combination 
with the measured ages of GCs in the Milky Way that 
instead constrain best their lower redshift formation rate, 
allow us to conclude that there are two distinct epochs 
of proto-GC formation: one at z ~ 2.5 near the peak of 
the global SFR density in the universe, and one before 
or around the reionization epoch. Figure [6] (left) shows 
the upper limits on the GCs formation rate (top panel), 
the proto-GC redshift distribution (middle panel) and the 
cumulative proto-GC redshift distribution (bottom panel) 



from the combined constraints from the LF and colors dis- 
cussed in the previous section. In addition, the proto-GC 
redshift distribution and cumulative distributions are con- 
strainecQ at re dshifts z < 3 by th e meas ured ages of Milky 
Way GCs from iForbes fc Bridges! (|201Ch . The lines in Fig- 
ure [5] (left) correspond to different assumptions on the min- 
imum mass of the halos in which proto-GC systems form: 
T mr > 5 x 10 4 K (dashed line), T vir > 8 x 10 4 K (dotted 
line); and T vir > 1.5 x 10 5 K (solid line). The limits on the 
GCs formation rate are independent of the assumed fiducial 
value for p gc , while the redshift distribution and cumulative 
redshift distribution of GC are normalized to p gc (that is 
somewhat uncertain as it depends on /<«)■ Because of the 
uncertainty on fdi and because the lines represent upper lim- 
its, the cumulative GC distribution may not be normalized 
to unity. 

Although our results only set upper limits on the GC 
formation rate, at redshifts z £ 6 the LFs and colors of sim- 
ulated proto-GC systems are in good agreement with ob- 
servations, suggesting that a large fraction of high redshift 
galaxies have a star formation mode compatible with proto- 
GC formation. Moreover, the upper limit on the cumula- 
tive distribution of proto-GCs normalized to p gc approaches 
unity, also suggesting that the upper limits are close to the 
actual values of the distribution. We also note that even a 
small fraction (10-20%) of present-day GCs forming at z £ 6 
accounts for a significant fraction of the total star formation 

1 This constrain implicitly assumes that Milky Way GCs accu- 
rately describe the age distribution of the GC population in the 
rest of the Universe. 
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activity at high redshift. Our model shows that if proto-GCs 
are allowed to form in halos with T V i T Si 5 x 10 4 K (dashed 
lines), current deep fields are missing a significant fraction 
of the global SFR because of the steep faint-end slope of the 
LF. 

Finally, we re-asses the possibility that proto-GCs ma y 
be dominant sources of IGM reionization (Ricotti 2002). 
The high star formation efficiency of young proto-GCs and 
their location in the outer parts of galaxies both point to 
a high escape fraction of the ionizing radiation they emit. 
Figure [6] (right) shows the upper limits on the number N 
of ionizing photons emitted per Hubble time per baryon by 
proto-GC formation, assuming f e3C — 1. The meaning of 
the lines is the same as in the left panel. A value of J\f ~ 5 
is sufficient for reionizing th e fGM and keeping it ionized 
l|Miralda-Escude et al.ll2000h . This scenario is especially in- 
teresting in light of current studies that are trying to identify 
the sources of reionization in deep surveys at z ~ 7 — 8, find- 
ing that the sources of reionization must be faint and thus 
relatively numerous dBouwens et aljfeoich . 



6 SUMMARY AND CONCLUSIONS 

Present day GCs are difficult to observe beyond the Milky 
Way because their stellar populations are very old and faint. 
All their bright massive stars have evolved beyond the main 
sequence and the surviving clusters have lost a significant 
fraction of their lower mass stars due to dynamical effects 
(tides, two-body relaxation, dynamical friction). However, 
during the first 5-fO Myr after their formation, proto-GCs 
were extremely bright, with the most massive GCs reaching 
absolute magnitudes of Mv ~ —19. This is due to their 
high- efficiency of star formation and their relatively large 
total mass at formation that can been estimated about 5- 
10 times larger than their current stellar mass [although 
this number is uncertain , depending on dynamical effects 

Ostriker fc Gncdin 1997) as well as the IMF of their stars 

Schaerer fc Charbonnelll201ll )]. 



In this paper we have shown that the most massive 
proto-GCs are detectable in the HUDF and will be de- 
tectable by the JWST up to redshift z ~ 8, when imaged at 
their peak luminosity. In addition, since GCs are found with 
a high specific frequency in present day dwarf galaxies (5-10 
within of a few kpc from the dwarf center), and we expect 
their formation to be synchronized to within ~ 10 Myrs, 
such proto-GC systems appear a poorly resolved high red- 
shift galaxy with multiple knots of star formation. The short 
dynamical time scales in primordial dwarf galaxies support 
this scenario. However, because proto-GCs maintain their 
peak luminosity for a very small fraction of the Hubble time, 
their numbers are small within the limited field of view of 
deep HST and JWST images, especially if a large fraction 
of them form at low redshifts when the Hubble time is large. 
But if a non- negligible fraction of today's GCs have formed 
at redshifts z > 3, they are detectable in the HUDF and 
will contribute a substantial fraction of the high redshift 
ste llar populatio ns and to IGM reionization, as first noted 
bv lRicottil i|2002l) . 

Our present work has shown that only a small fraction 
of the present-day GCs populations could have formed in the 
redshift range 3 < z < 5, because otherwise they would af- 



fect significantly and be inconsistent with the observed LFs 
and the colors at 1500 A of high-z galaxies in the HUDF. In 
addition, given the age estimates of Milky Way GCs, a sig- 
nificant old population that formed at z > 3 must exist (even 
if we take into account the large errors on their age estimate 
from color- magnitude methods). We conclude that the bulk 
of the old GC population with estimated ages ^ 12 Gyr, 
roughly 50% of the total, formed in the most massive dwarf 
galaxies at redshifts z ^ 6 and that proto-GC formation, an 
important mode of star formation in some primordial dwarf 
galaxies, likely dominated the reionization process. 

Although our study only sets upper limits on the for- 
mation rate of proto-GCs, the steep faint-end slope of the 
LFs and the blue colors of z ~ 7 — 8 galaxy candidates in 
the HUDF are consistent with a significant fraction of their 
stellar populations being proto-GCs. The simulated images 
of systems of proto-GCs forming in z ~ 7 — 8 dwarf galax- 
ies present remarkable morphological similarities with some 
of the observed candidates. JWST will be able to detect 
and image this population that is likely responsible for the 
reionization of the universe, to higher redshifts and fainter 
magnitudes. 

Finally, we speculate that if GC formation was a dom- 
inant mode of star formation in some dwarf galaxies, this 
may have implications to understanding s ome unsolved is- 
sues r egarding the Milky Way satellites. iBovill fc Ricottil 
(2009) have pointed out that many of the ultra-faint Milky 
Way satellites have properties that are consistent with be- 
ing the fossil remnants of the first dwarf galaxies formed 
before reionization. But they also noted that a prob- 
lem remains regarding an apparent lack of massive and 
bright satellites that cannot be reconciled by simply pop- 
ulating the most massive sub-halos with luminous satel- 
lites uarig_jnonotonic^^ meth- 
ods (|Bovill fc RicottteOUal lbTl. lBovlan-Kolchin et all (|201ll ) 
have also suggested a possible problem with the density 
of the most massive Milky Way satellites that seems in- 
consistent with dark matter ACDM simulations (neglect- 
ing baryonic physics). Solutions to this second problem in- 
clude SN feedback processes that may reduce the central 
density of satellites, and statistical arguments involving a 
better knowledge of the mass and circular velocity of the 
Milky Way. However, thes e two solutions would no t solve 
the problem emerging from IBovill fc Ricottil l|2011bl ) work. 

A solution that would reconcile both results is to as- 
sume that classical dwarfs do not reside in the most mas- 
sive dark matter satellites. This is possible only if the most 
massive satellites of the Milky Way are darker than some of 
the least massive (i.e., a non-monotonic mass-to-light ratio). 
There are two ways to achieve this result: i) suppress star 
formation preferentially in the most massive satellite pro- 
genitors, with a feedback process such as photoheating by 
AGN etc.; or ii) remove most of its stars without removing 
most of the dark matter. This second proposal seems im- 
plausible if stars are concentrated at the center of the dark 
matter halo. However, if star formation in some of the most 
massive surviving Milky Way satellites was dominated by 
proto-GC formation in the outskirts of their dark matter 
halo (as observed in nearby dwarf galaxies), tidal stripping 
of their GC systems may significantly reduce their luminos- 
ity without destroying the dark halo, thus explaining their 
apparent darkness. Examples of galaxies with stellar popu- 
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Figure 6. Left. Upper limits on the GC formation rate as a function of redshift from LF and color constraints for systems of GCs 
(see text). The top panel shows the GC formation rate compared to the observed SFR (long-dashed line), the middle panel shows the 
distribution of formation redshifts of GCs and the bottom panel shows the cumulative distribution. The different lines correspond to 
different assumptions on the minimum virial temperature the host halos in which GC systems are allowed to form: T„j r > 5 X 10 4 K 
(dashed line); T v i T > 8 X 10 4 K (dotted line) T v i r > 1.5 X 10 5 K (solid line). The thick lines include the constraints at low-redshift from 
estimated GCs ages. Right. Upper limits on the ionization rate per Hubble time per baryon from GCs. The lines correspond to the GCs 
formation rates as in the left panel. 



lations clearly dominated by GCs have not been identified 
in isolation in the local universe, but their formation might 
be biased, such that today they have mostly merged into 
massive galaxies. Another possibility is that isolated mas- 
sive dwarfs were once dominated by GC systems b ut have 
accre ted gas and formed new stars at later times (|Ricottil 
2009), while the ones that merged into the Milky Way did 
not experience further star formation. 
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